COMBINATORIAL TOPOLOGY OF THREE-DIMENSIONAL 

SELF-AFFINE TILES 



CHRISTOPH BANDT 

Abstract. We develop tools to study the topology and geometry of self-affine 
fractals in dimension three and higher. We use the self-affine structure and obtain 
rather detailed information about the connectedness of interior and boundary 
sets, and on the dimensions and intersections of boundary sets. As an application, 
we describe in algebraic terms the polyhedral structure of the six fractal three- 
dimensional twindragons. Only two of them can be homeomorphic to a ball but 
even these have faces which are not homeomorphic to a disk. 



1. Introduction 

When Levy [32J introduced his famous curve in 1938, he also constructed fractal 
surfaces in a similar way. 70 years later, we have plenty of papers on fractal sets 
in the plane, and a number of general statements and constructions which hold in 
every dimension, but very few studies on the geometry of fractals in dimension > 3. 

Visualization is certainly one reason to prefer plane sets: they can be easily shown 
on a computer screen. In dimension 3, visualization is more difficult, even though 
there are good ray-tracing programs, like chaoscope [15] and IFS builder [27] , which 
allow to look at three-dimensional fractals from any chosen viewpoint, with pre- 
scribed light sources. Figure [1] shows two views of one of our main examples. Does 
it have interior points? If so, is the interior connected? Is the figure simply con- 
nected? Unlike in two dimensions, the answer to such questions can hardly be 
guessed from the pictures - this requires tools which we develop here. 

There are several mathematical difficulties in dimension > 3. We avoid the prob- 
lem that linear maps do usually not commute, by concentrating ourselves to fractals 
generated by a single matrix. But there is another problem which will concern us: 
proper self-similarity is rather an exceptional case. The moduli of all eigenvalues of 
a 3 x 3 matrix do rarely coincide while for a 2 x 2 matrix they coincide as soon as 
the eigenvalues are complex. This problem is addressed in Section 2, and Theorem 
2.3 implies that in M 3 there is essentially only one type of 'integer' self-similar tile 
with 2 up to 7 pieces. 

Other difficulties concern the topology of three-dimensional fractals. Our main 
aim here is to study the topology of self-affine tiles, and of their boundaries which 
are fractal surfaces. In particular, we would like to know which of these tiles are 
homeomorphic to a closed ball. In two dimensions, the corresponding question for 
disk-like tiles was discussed in a large number of papers HI [FJ [121 EES E21 EE EH EH] 
(more references can be found in [311 021 S3]), and Jordan curve arguments were 
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Figure 1. Two views of the twindragon C , made by chaoscope [15j. 
Is this set homeomorphic to a ball? Is the interior connected? We 
develop concepts and tools to answer such questions. 

heavily used. Among others, it is easy to see that connectedness of the interior of 
a self-affine tile in the plane is sufficient for its disk-likeness. In M 3 this is not the 
case. The topology of the boundary of a tile can be much more complicated than 
in M 2 . Even the connectedness of a tile is not so easy to verify |JJ [26] . 

The geometric phenomena which we encounter in fractal tiles are not caused by 
knots or wild spheres [H]. It is rather the self-affine fibre structure of the boundary 
which makes the topology complicated. Apparently, the different eigenvalues of 
the generating matrix produce long and thin fibres which can pierce the interior of 
neighboring tiles, or distort the boundary structure. 

We develop algebraic tools which describe the geometry of a self-affine tile T in 
arbitrary dimension, in a similar way as homotopy and homology groups describe 
the geometry of manifolds. Our tools use self-affinity in a specific way, and they yield 
quite detailed information about the geometry. The basic concept is the neighbor 
graph G = (V,E). It can be considered as a blueprint of T which contains all 
information about the topology of T in a simple scheme. In the case where T 
tiles by translations in a unique way, the vertices v G V are those translations for 
which T + v belongs to the tiling and intersects T. Such translations v will be called 
neighbor maps, and they do not only represent the neighboring piece T + v but also 
the boundary set T D (T + v) of T. This set can be a face, an 'edge' or just a point 
of T, or a more complicated fractal boundary set. For a ball-like tile we expect that 
faces are homeomorphic to disks and edges are homeomorphic to intervals, and we 
have to check this idea. 

Neighbors in lattice tilings were introduced by Indlekofer, Katai and Racsko [23] , 
and neighbor maps for self-similar sets were defined by Bandt and Graf |7J. For a 
self-affine tile, the boundary sets have themselves a self-affine structure: each face 
or edge is the union of smaller copies of other boundary sets. This is indicated by 
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the edges of the graph G. The definition of G in Section 3 implies that the boundary 
dT is a self-afline graph-directed construction in the sense of Mauldin and Williams 
[34j, and the graph which directs this construction is G itself. This fact was stated 
by Scheicher and Thuswaldner [36] for lattice tiles, and used in many other papers 
to determine the Hausdorff dimension and topological structure of dT for plane self- 
similar tiles [211 E3 [251 |33l EB EH |39]. Some authors [161 EE El HJ [2] prefer the 
contact matrix defined by Grochenig and Haas [22] which describes a subgraph of 
G. 

A first application of the neighbor graph concerns the existence of self-affine lattice 
tilings, a problem treated in [22J and in a series of papers by Vince [101 El 112] . 
Theorem 5.1 shows that the existence of such tilings is equivalent to the compatibility 
of all neighbors, which can be easily checked with G. 

In sections 7 to 12 we study the topology of the boundary of T. In section 8 we 
define the hierarchy of different boundary sets like faces, edges, points etc. which 
corresponds to the hierarchy of sets with different Hausdorff dimensions in Mauldin- 
Williams constructions. An important problem is to determine all faces - those 
boundary sets which cover an open set on the boundary. Two methods are developed 
for this purpose. Theorem 9.2 gives an algorithmic approach which only uses the 
combinatorics of the neighbor graph. Theorem 10.1 provides an analytic approach 
based on recent work by He and Lau [23] and Akiyama and Loridant [2 J , connected 
with eigenvalues of adjacency matrices and Hausdorff dimension. In section 11 we 
define a polyhedral structure of fractal tiles by the intersection sets of two and more 
faces. There are corresponding neighbor intersection graphs G 2 ,G 3 ,... which allow 
to calculate details of the polyhedral structure and compare it with the structure of 
ordinary polyhedra. 

In contrast to the contact matrix, all these tools can also be applied to tiles with 
incompatible neighbors as well as to fractals which are not tiles, and neighbors can 
be related by arbitrary isometries of 1R 3 instead of translations. To illustrate our 
methods, we selected a small class of examples, which is introduced in section 6: 
the twindragons. These are the self-affine lattice tiles with two pieces. The one- 
dimensional twindragon is the interval. In the plane, there are three examples: 
rectangle or parallelogram, the twindragon and the tame twindragon. They are 
all disk-like [2H HJ- In three-dimensional space, there are seven examples, which 
we denote by letters A to Q . While A is conjugate to a cube, T and Q have an 
extremely intricate structure. We mention some of their properties without proof. 
For T> and £ we determine the boundary faces and their intersections, and we show 
that their interior is not simply connected. B and C seem to be homeomorphic to 
a ball. In section 12 we determine their exact polyhedral structure. We show that 
they have faces which are not homeomorphic to a disk, but C has connected interior. 

This paper focusses on concepts rather than algorithmic questions although it 
seems clear that the neighbor graph and related finite automata are well-suited for 
computerized evaluation. Section [TJ] gives a short account of algorithmic aspects. 
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2. Self-affine LATTICE TILES 

Basic concepts. The following definitions are standard [28j ESI H21 03]. A 

lattice L in M. n is the set of integer linear combinations of n linearly independent 
vectors ei, ...,e n . Usually we take L = Z n . A regular-closed set T is called a tile if 
it admits a tiling of lR n . A tiling by the tile T is a countable union Ufc^fc which 
covers R™, such that each is isometric to T, and no two tiles have interior points 
in common. We have a lattice tiling by L if the tiles Tjt are just the translates T + k 
with G L. 

A tile T is a self-affine lattice tile if there is an affine expanding mapping g and 
a lattice L such that g preserves L and maps T to a union of tiles T + ki with 
ki G L. With respect to the basis vectors e«, the map g has a matrix representation 
g(x) = Mx where M is an integer matrix. Expanding means that all eigenvalues 
of M have modulus greater one. T is determined by the map g and by the lattice 
coordinates ki, k m of the tiles which form the supertile g(T). 

m 

g(T) = MT={jT + k j (1) 

i=i 

For T to become a tile, we must have m = | detM|, and the set K = {ki, k m } 
must fulfil some condition. A sufficient condition is [5] 

g(L) = {jL + k j = L + K (2) 

3=1 

in which case K is called a standard digit set or complete residue system. The latter 
name comes from the fact that for | detg| = m, the subgroup g{L) of the additive 
group L has m residue classes, so by ([2]) each class is represented by exactly one ki. 
See Lagarias and Wang EH] who also studied non-standard digit sets. In this 
paper, we consider only standard digit sets. 

The self-affine tile T is called s elf- similar if the mapping g is a similarity mapping 
with respect to the Euclidean metric. In this case, we must distinguish the standard 
basis which defines the Euclidean distance, and the basis for the lattice L which 
need not be orthonormal. 

Conjugacy of tiles. T is said to be conjugate to a self-similar tile if there is 
a linear map h so that T = h(T) is a self-similar tile. The expanding map for T is 
g = hgh' 1 , the lattice is L = h(L), and the digits are ki = h(ki), for which ([T]) and 
([2]) are easily checked. 

Example 2.1. In the plane, we consider the digits k\ = (°), &2 = (q)- The parallelo- 
gram T with vertices (°), ( ), ( x ), (^) is a self-similar tile with respect to M = (^J 
and the lattice L = 1? . The rectangle T' = [0, 1] x [0, a/2] is also a self-similar lattice 
tile, with M' = and L' = {(J^^ \ nt, n G Z}. The unit square T" = [0, l] 2 is 

a self-affine tile with respect to the matrix M" = (° 2 ) and L" = Z 2 , but not self- 
similar. However, T" is conjugate to the self-similar tile T' since h(x) = (q 1 ^)^ 
fulfils h(T") = T'. Note that h{T) = V for h(x) = 
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We shall identify conjugate tiles so that T,T',T" are considered as different ver- 
sions of the same tile, which is essentially self-similar. Considering tiles up to con- 
jugacy means that we focus on the essential data: the characteristic polynomial of g 
instead of different matrix representations. We shall also assume that h = (J) since 
this can be obtained by conjugacy with a translation. 

Proposition 2.2. (The class of self -similar tiles) 

A self-affine lattice tile T with respect to g{x) = Mx is conjugate to a self-similar 
tile if and only if all eigenvalues of M have the same modulus. 

Proof. Necessity of the condition follows from the fact that eigenvalues are not 
changed under conjugacy. On the other hand, if the eigenvalues have equal modulus, 
there is a matrix B such that g(x) = BMB~ 1 x is a similarity map. Thus h(x) = Bx 
maps T to a self-similar tile. □ 

Theorem 2.3. (Very few self-similar lattice tiles in dimension 3) 

If in M. 3 a self-affine lattice tile T with respect to g(x) = Mx is conjugate to a 

self-similar tile, then one of the following two conditions is true. 

(i) m — | det M\ is a cubic number - in particular m > 8, 

/ ±m 

(ii) M is conjugate to M — I 1 

\ 1 

Proof. With respect to the lattice base, g has an integer matrix, and so the 
characteristic polynomial p(\) of M has integer coefficients. We express them in 
terms of the eigenvalues Aj. 

p(X) = -A 3 + (Ai + A 2 + A 3 ) • A 2 - (AiA 2 + AiA 3 + A 2 A 3 ) • A + AiA 2 A 3 

where AiA 2 A 3 = det M = ±m. Since T is conjugate to a self-similar tile, all Aj have 
equal modulus: |Aj| = r > 1. 

If the eigenvalues are real, then A« = ±r and m = ±r 3 . In this case Ai + A 2 + A 3 G 
{±r, ±2r, ±3r}, so r must be an integer and m a cubic number. 

Now let us assume Ai, A 2 are complex eigenvalues. Then Ai = r(cosa; + isina) 
and A 2 = r(coso; — isina) for some a, so 

p(A) = -A 3 + 2r (cos a + |) • A 2 - (r 2 + 2sr 2 cos a) ■ A + r 3 s 

where s = ±1 is the sign of A 3 . The coefficient of A can be written as — 2sr 2 (cos 
Since the coefficients are integers, either r must be a rational number - and hence 
an integer, and m a cubic number. Or cos a + | must be zero. In this case, p(A) = 
—A 3 + sin which is the characteristic function of M. □ 

Remark 2.4. The orthogonal part of M is a rotation around 120° for s = +1, and 
a rotation around 60° composed with a reflection at the plane of rotation for s = —1. 
Even in the case where r is an integer and m a cubic number, 2r cos a must be an 
integer and a can assume only few values. Thus there are really very few self-similar 
lattice tiles in M, 3 . 
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3. Neighbors of tiles and self-affine sets 



When we want to study the boundary B of a lattice tile T, it is quite natural to 
consider the neighbors in the tiling - those tiles T + k for which = T fl (T + k) 
is non-empty. The B^ can be considered as the faces of T, and B is the union of 
the B^. Obviously, the number of neighbors is finite. This idea was introduced by 
Gilbert [21] and Indlekofer, Katai and Racsko [21] and used in many other papers. 

Neighbors in fractals. There is a related concept of neighbor in self-similar 
fractals which appeared first in [7]. Suppose we have a fractal A which consists of two 
copies or itself: A = f\{A) U f 2 {A), where the fi denote contracting similarity maps. 
Then the geometry and topology is determined by the structure of the intersection 
set C = fi(A) fl /2(A). Only at points of C it makes sense to zoom into the picture. 
If x is in fi(A), say, and U is a neighborhood of x in A which does not intersect C, 
then the magnified copy fx (U) of U does already exist in A. 

The sets /{" (C) and f 2 1 (C) can be considered as boundary sets of A, where 
A intersects a potential neighbor f^ 1 f 2 (A) or f 2 fi(A). These need not be the 
only boundary sets, however, since C consists of intersections of smaller pieces like 
/1/2 (-4) 0/2/2(^4), and these may be at other positions when they are pulled back to 
A. This argument shows that the boundary sets have a self-similar structure. The 
number of boundary sets obtained by going to smaller and smaller pieces need not 
be finite, but in most common examples it is. 

Let us define neighbors in this general sense. Let fi, f m denote contractive 
affine mappings on R n . There is a unique compact non-empty set A with 



the self-affine set with respect to the fi, cf. [T7]. For each integer q the set A splits 
into m q small copies fj(A) = Aj, where j = ji-..j q denotes a word of length q from 
the alphabet I = {1, m}, and /j = fj 1 ■ ... ■ fj . 

A potential neighbor of A has the form h(A) = f^fj{A) where i,j G I* are any 
words over I such that the pieces A\ = fi(A) and Aj = fj(A) intersect. The idea 
is that f-' 1 maps A{ to A and Aj to h(A), a neighbor set which intersects A and 
should be of comparable size, h is called a neighbor map and the set B = A fl h(A) 
the corresponding boundary set of A. To get 'comparable size', we confine ourselves 
to words i, j of equal length and make some assumptions concerning the fi. 

Proposition 3.1. (Neighbor maps which are isometries) 

Let A C M. n be a self-affine set of the form Let i, j be words of the same length 
q, and h = fr l f- y 

(i) If all fj are similarity maps with the same similarity factor r, then h is a 
Euclidean isometry. 

(ii) If fj(x) = M~ 1 (x + kj) where M is an expanding matrix and ki,...,k m G M n , 
then h is a translation. 



rn 




(3) 
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Proof, (i): If the fj are similarity maps with factor r then /j and f^ 1 are similarity 
maps with factor r q and r~ q , respectively. Thus h is an isometry. 
(ii) is first proved for q = 1. Since f^ 1 {y) = My — ki, we have 

fr l f j (x)=x + k j -k i (4) 

Moreover, if h(x) = x + v is a translation then 

fr 1 hf j (x)=x + Mv + k j -k i (5) 

Now induction on g shows that for /i = /j - /j, 

= /, /•••/,: ■// •-./),(•'•) = x + M^fe - fcij + M"" 2 (% 2 - fc l2 ) + ... + - fcj 

(6) 
□ 

Lattice tiles as a special case. The case (ii) includes our self-affine lattice 
tiles since the defining equation (jTj) can be rewritten as 

m m 

T = {J g-\T) + g~ l {kj) = |J M-!(T + fc,-) (7) 

3=1 3=1 

We get fj(x) = M _1 (x + fc^) as contracting maps for the self-affine set T. 

At this point it is possible to explain the concept of standard digit set ki,...,k m 
which says that kj — k iq is not in ML for j q ^ i q (see section 2). This property 
implies that the translation vector in can never be zero if j q ^ i q since kj — ki 
cannot cancel with the other terms in (jSJ). As a consequence, the inequality h ^ id 
then holds whenever i ^ j. From the results in \7\ then follows the so-called open 
set condition which says that the pieces f\{T) and fj(T) have no interior points in 
common. (In [5], this was shown by Baire's category theorem.) In other words, T 
and a translate h(T) in ([6]) have no common interior points. 

4. The neighbor graph 

Self-similarity of boundary sets. To explain the self-similar structure of 
boundary sets, we return to the general case. Consider intersecting pieces Ai and 
A- } in the self-affine set A. Since all pieces divide into m subpieces, like A in (J3]), 

rn 

AnAj= |J Au n A$j. 

i,3=l 

Now consider the boundary sets B = /-^(Ai D Ay) and = /-^(Au fl A- S j). Then 

tm \ m m 

|J AunAjA = |J fifr'fr^AunAu) = [j / 4 (By). 
ij'=l / t,j=l i,j=l 

At B,Bij we suppressed the subscripts i,j since such an equation holds for each 
possible boundary set. In other words, the subdivision of pieces induces self-similar 
representations of the type d3J) for all boundary sets - not with one type of set, but 
with several types. We must assume now that we have only finitely many possible 
boundary sets, which is true for lattice tiles with standard digit sets. The unions 



8 



CHRISTOPH BANDT 



on the right contain a lot of empty terms, so we introduce a graph which better 
describes the system of equations. 

Concept and properties of neighbor graph. The neighbor graph G = (V, E) 
of a self-affine fractal A has as vertex set all neighbor maps h = f^ l f- 3 , and a directed 
edge marked with i goes from htoh' — f^ 1 f^ 1 fjfj- Loops and multiple edges (even 
with the same label) are possible [10]. For a formal definition, we focus on self-affine 
tiles. 

Definition 4.1. LetT = [J fj(T) be a self-affine lattice tile with fj(x) = M~ l {x J rkj) 
so that the neighbor maps have the form h(x) = x + k. We identify h with the 
translation vector k, and asopciate it with the boundary set Bk = T R (T + k). An 
edge from k to k' with label i is drawn when relation |3J) holds: 

G = (V, E) with V = {k\Tf](T + k) ^ 0} and 
in 

E = [jE l with Ei = {(k, k', i)\k' = Mk + % - k,} (8) 
i=i 

Here Ei denotes the set of edges with label i. It should be mentioned that V 
contains a root vertex k = (or h = id in the more general notation) which does 
not correspond to a boundary set. There are edges (0, kj — ki, i) from the root which 
correspond to the first step (jl]) in the calculation of neighbor maps. The loops 
(0, 0, i) will not be drawn, however. For the calculation of G, see section [TU [101 [11] 
and the example below. 

Proposition 4.2. (The boundary equations) 

Let T = Ujli fj{T) be a self-affine lattice tile with neighbor graph G = (V, E). The 
boundary sets Bk corresponding to the k G V \ {0} fulfil the following equations. 

B k = \J UA B k>\j e {l,-,m}, (k,k'J) G E} (9) 

This is an immediate consequence of Definition 14.11 and of the discussion above. 
Mauldin and Williams [31] called families of such fractals Bk graph- directed con- 
structions, and proved that the Bk are uniquely determined by the graph G and the 
maps fj. At the end of section 10 we briefly discuss the open set condition of (J9]). 

Proposition 4.3. (The basic symmetry of G) 

If h is a neighbor map, then h~ l is also a neighbor map. For translation this means 
that to every k there is a —k, and each boundary set Bk = T D (T + k) has an 
opposite set 

B_ k = {T - fc) n T = B k -k. 
This symmetry of V will extend as a graph isomorphism to the edges, but not to the 
labels. If k' = Mk + kj — ki then (k, k') has label i and (—k, —k') has label j. 

In a previous publication with M. Mesing [10], we wrote both labels i,j at the 
edge (k,k'). This is not necessary here since we can always recover the second label 
from the opposite edge. 

A two-dimensional example. In our examples, neighbors will be denoted by 
lower-case Roman letters, and — b will denote the opposite vertex of b. We start 
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Figure 2. A two-dimensional tile and its neighbor translations. 

with a two-dimensional tile [5j H2] which can be considered as a modification of the 
square, and as an extension of the Sierpihski gasket. The construction of neighbor 
graphs for three-dimensional tiles proceeds in the same way as demonstrated here. 

Example 4.4. Let fj(x) = (x + kj)/2 for j = 1, ...,4 where k\ = Q, k 2 = = 
(°), and ki = (ij. The tile T is shown in Figure® Clearly, T C [—1, +1] 2 . Thus a 
translation x = (^) with T(l(T + x) ^ must fulfil \x±\ < 2 and \x 2 \ < 2. Moreover, 
L is the integer lattice, so x\,x 2 G {—2,-1,0,1,2}. In this way we prove that for 
an arbitrary self-affine lattice tile, the set V is finite. 

The tile T has six neighbors which intersect T in a single point, and six neighbors 
which intersect T in an uncountable set which is in fact a Sierpinski gasket. 

Proof. We start with the translation vectors k = kj — k,i where i ^ j, see the 
right-hand part of Figure [2j If k is in V, there is the edge (0, k,i). Since our matrix 
is M = (q 2)5 a ll these points do belong to V, with a loop (k, k,j) at each point k, 
because kj k^ = *2(^kj k^ ~\~ ki kj. 

Moreover, these will be the only points of V. For the other points x with < 
2, \x 2 1 < 2 the recursion ([5]) will lead to vectors k' which are outside the range of 
Xi,x 2 and thus are no neighbor maps. A simple calculation [TO] shows that once we 
are outside the range, we can never come in again by the recursion ([5]). 

Thus the translations ±a = ±(^),±b = ±(" 1 ),±c = ±(_° 1 ) marked by •, and 
±( 2 ), ±(2), if^i) marked by o are the vertices of our neighbor graph. The last six 
vertices have only the loop and no further outgoing edges. These six vectors are the 
translations along the sides of the big triangle in Figure [21 They translate T to a 
neighbor h(T) which has a single point as intersection with T. These neighbors are 
called point neighbors. 

To find point neighbors, we need only look at the graph, not at the picture. The 
loop at k — ( 2 ) with label 2 says that the boundary set has address 2222... = 2, 
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Figure 3. The reduced neighbor graph for the Sierpihski tile. 

and the opposite vector has address 4. This means that T fl (T + k) is exactly 
one point, which has address 2 in T and address 4 in T + k. Indeed the intersection 
is the point k 2 = ( ) which is the fixed point of f 2) and T + k will meet this point 
with the translate of k±, the fixed point of f^. The translate of k 3 along the vector 
k' = (~ ) is also k 2 . This way it turns out that two of the six point neighbors will 
meet T at each vertex of the triangle. 

We now forget the point neighbors o. The graph of the remaining six • neighbors 
±a, ±b, ±c and the root, shown in Figure EJ will be called the reduced neighbor 
graph [I0]. On the left of Figure EJ the corresponding boundary sets are indicated. 
As mentioned above, each vertex k = kj — ki has a loop with label i. The edge 
(— b, a) has label 2 and the opposite edge has label 3 because a = 2 ■ (— b) + /c 3 — k 2 . 
In a similar way, all other edges were determined. 

The graph in Figure [2] is irreducible which means that all corresponding boundary 
sets have a similar structure. This also follows from the equations of the graph- 
directed system ([9]): 

#a = /i(B a )U/ 2 (B_ c )U/ 3 (B_ b ), B_ h = f 2 (B a UB_ h UB c ), £L C = f 3 (B a UB h UB_ c ). 

All these boundary sets are in fact small Sierpihski gaskets, as can be seen in the 
figure. The theorem below shows that they are uncountable. □ 

Topology and address map. A sequence B = Sis 2 ... G I°° is called address of 
the point 1 in a self-afline tile T if it is contained in the pieces T n ___ Sn for n — 1,2, ... 
(cf. [13J, Chapter IV, [T7], Section 9.1). The points of the intersection sets Tj fl Tj 
have two addresses, one starting with i and the other one with j. The address map 
71 : I°° ->Tisa quotient map, it determines the topology of T. For that reason it is 
important to know which addresses will be identified. The neighbor graph provides 
this information, as can be easily shown ([10], sections 4 and 5). A path in the 
directed graph G is a finite or infinite sequence eie2... such that the terminal vertex 
of coincides with the initial vertex of e^+i. A finite path cycle if the 

terminal vertex of e n is the initial vertex of e\. 
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Theorem 4.5. (Neighbor graph and addresses of points, see [lOj ) 

Let T be a self-affine tile with the address map it : I°° — > T and the neighbor graph 

G=(V,E). 

(i) Two sequences 6 = siS2-.. and t = tit 2 ... are addresses of the same point if 
and only if there is a path e^... starting in the root with labels siS2---, such 
that the opposite path -e\-ei-.. has labels t\t%... 

(ii) The addresses of the points of a boundary set coincide with label sequences 
of paths starting in the vertex k. They form a regular language 

(iii) If only one cycle can be reached from the vertex k by a directed path, then 
By. is a singleton. 

(iv) If two different cycles can be reached from k, and these cycles can be reached 
from each other, then the set B^ is uncountable. 

Note that a cycle with a diagonal, or with a loop at one of its points counts as 
two different cycles which can reach each other. 

5. Existence of self-affine lattice tilings 

Before we can go on, we must settle one difficulty with neighbors. They need not 
be compatible: h(T) and h(T) can be neighbors of T, but not of each other. In that 
case, different tiles in a tiling will have different neighborhoods. 

It turns out that the compatibility of neighbors is connected with another question 
which was studied by Grochenig and Haas [22] and in several papers of Vince [40j 
HU H2] as a "central" problem in the field. A self-affine lattice tile with standard 
digit set will always admit self-affine tilings. It is a non-trivial result that it will 
also admit lattice tilings [SU] but these need not be self-affine in the sense that tiles 
assemble to form supertiles ([12], Example 4.2). Vince ([2], [12], Theorem 4.3) gave 
10 equivalent conditions for the existence of a self-affine lattice tiling. Two of them 
are of an algorithmic nature. 

Here we extend Vince's list by proving that a self-affine lattice tiling exists if and 
only if all neighbors are compatible with each other, and we show that this property 
can be easily decided with the neighbor graph. The following theorem also shows 
that all self-affine sets with mappings fj{x) = M~ l (x + kj) which are not lattice 
tiles must have incompatible neighbors. 

Theorem 5.1. (Compatible neighbors and self-affine lattice tilings) 
Let A C M n be a connected self-affine set with respect to mappings fj{x) = M~ 1 (x + 
kj),j = 1, ...,m where M is an expanding matrix and the open set condition holds. 
Then the following conditions are equivalent. 

(i) All potential neighbors appear together at one piece A\. 

(ii) A is a self-affine lattice tile which admits a self-affine tiling by a lattice. 

(iii) In the neighbor graph, each vertex k ^ has at least one incoming edge with 
each of the labels j = 1, m. 

Proof. (i)=^(ii): We assume A is a self-affine set, and A± has all possible neighbors. 
Then the number of neighbors must be finite, at most m n if i = i\...i n . Moreover, 
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consists of copies of A which include all possible neighbors h(A). Since the 
maps are translations h(x) = x + k by Proposition 3.1, let again K denote the set 
of all these translation vectors. We define an infinite pattern of translates of A as a 
union of increasing compact sets: 

oo 
9=1 

All translates A' of A in this pattern have the form f^ q {A\- } ) for some number q 
and some word j. This implies that they all have the same neighbors as A. To 
see this, consider the (potential) boundary of a subpiece A^ with j G {l,...,m}, 
defined as the union of all its intersections with potential neighbors. This boundary 
is contained in the union of the boundary of Ai and the intersections A-^ D Ay with 
i — 1, m, i 7^ j. So by assumption the boundary of A-^ is contained in A, and by 
induction on the length of j this is proved for subpieces Ay. 

Thus all 'atoms' A 1 of our infinite pattern have a maximal set of neighbors which 
covers their boundary completely. But there is only one maximal set: the complete 
set of translates A' + k with k G K. 

Let L denote the lattice generated by K. Since A was connected, any two pieces 
of the same level are connected by a chain of neighboring pieces (cf. [TTj, 8.2.1), and 
each atom has the form A' = A + k with fcei. Moreover, L has the same rank as 
the linear subspace generated by A, which we will now assume is M n . 

Since all A' have the neighbor vectors K, the infinite pattern coincides with A + L. 
As a consequence, the pattern must be a tiling of M n : If there was a small open set U 
outside A + L then there would be arbitrary large copies f[ q {U) outside A + L which 
would contradict the finite mesh size of L. Thus we have constructed a self-affine 
lattice tiling by copies of A. 

(ii) =^(iii): Let W 1 = A + L be a self-affine lattice tiling by the tile A. Then each 
tile of the tiling has the same set of neighbor translations K. Thus if k is in K, 
then k appears as a neighbor translation of any piece Tj of a supertile T, where 
j G {1, ...,m}. If the neighbor is another piece T of T, there is the edge (0, k,j) in 
G. If the neighbor of Tj is in another supertile T' then we have the edge (k',k,j) 
where k' is the neighbor map between T and T'. Thus for each vertex k of the vertex 
set K, there is an incoming edge with label j. 

(iii) =>(i): We assume that each vertex k in the neighbor graph G = (V,E) has 
incoming edges with each label j. We have to find a word i = such that Ai has 
all neighbors k G V. For each k there must be a suffix i n tk)--^N of i which consists 
of labels of a path from to k. 

It is not difficult to construct such a word i by induction, going paths backwards 
(against the direction of edges) towards 0. Start with a vertex k± and let ji = ji-.-j qi 
be the labels of a backward path from k± to 0. According to our assumption, there are 
also backward paths with label sequence ji from all other vertices. Their endpoints, 
which are different from 0, will form a set V\. 

Next, take a k 2 G V\ and a backward path j 2 = j qi +i---j q2 from A; 2 to 0. Also take 
backward paths from all other k G V\ with label sequence j 2 and denote the set of 
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their non-zero endpoints by V 2 . Continue with k 3 e V 2 , and so on. Since V n+ \ has less 
points than V n , we will have V n = for some n. The sequence i = j qn ...j qn _ 1 ...j qi ...ji 
will have all required suffixes. □ 

All examples considered in this paper will fulfil the condition of theorem 5.1. For 
Example 4.4 this can be seen in Figure [31 for the twindragons in the figures below. 



6. The seven three-dimensional twindragons 

Twindragons and their symmetry. We now define the family of examples 
we are going to study here. A twindragon is a self-affme lattice tile with m = 2 
pieces. The expanding matrix M has determinant ±2. In dimension 1, a twindragon 
is an interval. In IR 2 , there are three examples: the rectangle (Example 2.1), the 
ordinary twindragon and the tame twindragon, and they are all disk-like [7]. 

Since we identified conjugate tiles, we could take k\ = in section 2. For twin- 
dragons in W 1 , we can also choose k 2 arbitrarily, as long as k 2 is not in an invariant 
linear subspace of M. For the basis {k 2 , 9^2) , 5 ,T1_1 (^2)} the matrix M of g will 
always be the same, determined by the characteristic polynomial of g (see [20J and 
the proof of Theorem 6.2 below). For our twindragons in M 3 , we choose another 
normalization: 

fci = (-^,0,0)', ^ = (1,0,0)' (10) 

This has the effect that is the symmetry center of T, and in the neighbor graph, 
an edge with label 1 passes from the root to a = (1, 0, 0)'. 

Proposition 6.1. (Symmetry of twindragons) 

Each twindragon T in lR n has a symmetry center at c = \{k\ + k 2 ). The point 
reflection at c interchanges the two pieces Ti,T 2 , and each boundary set with the 
opposite set Moreover, each boundary set of T also has a symmetry center 
at c k = c+ |. 

Proof. We can take k\, k 2 from ( flOl) since the symmetry is not changed by an afline 
conjugacy. Thus f\(x) = M~ l x + k\, f 2 (x) = M~ x x — k\. With <p(x) = —x we obtain 
4>fi = h4> or fi = 0/3-i0 _1 since 0" 1 = 0. For any word i = i x ...i n on I = {1,2}, 
we have fi = where each j q = 3 — i q is the opposite symbol. Since the point 

of T with address fi = i\i 2 ... can be represented as x = lim^oo fi 1 fi 2 ---fi n {z) where 
z is any starting point [17J, this implies that transforms this point into the point 
with the opposite address. In particular 4>(T) = T and 0(Tj) = T^i from which it 
follows that cf) maps T\ fl T 2 onto itself. 

For each boundary set B^ = TC\(T+k), the map (f) transforms T + k to <fi(T) — k = 
T — k, thus (j)(Bk) = -B_fc. On the other hand we had B_ k — B^ — k, so ip(x) = —x + k 
maps Bk to itself, ip is the point reflection at c k = §. □ 

The seven twindragons. The last proposition shows that the tile generated 
by fi, f 2 can also be generated by fi<p and f 2 <f>, where is the point reflection at c. 
Since in M 3 the map has determinant -1, this has the consequence that we need 
only consider the matrices M with det M = 2. With -2 we get the same twindragons. 
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Figure 4. Twindragon B , shown first in [20] 

Theorem 6.2. (List of three-dimensional twindragons, cf. [26| [T] ) 

Up to conjugacy, there are exactly seven twindragons, each of which is uniquely 
determined by the pair of coefficients (a, b) of its characteristic polynomial 

p(A) = -A 3 + a\ 2 + bX + 2 : 



Twindragon 


A 


B 


C 


V 


£ 


T 


Q 


Parameter 


(0,0) 


(-1,1) 


(1,-1) 


(0,1) 


(2,-2) 


(1,0) 


(0,2) 



Proof. These are the only integers (a, b) for which p(X) has only roots of modulus 
> 1 HH Q]. For fci = 0, A; 2 = (1, 0, 0)' and the basis {k 2 , Mk 2 , M 2 k 2 } the affine map 

/ 2 \ 

g{x) = 1 b • X 

\ 1 a ) 

has characteristic polynomial p(X), and yields a twindragon. The standardization 
ffTUj) translates this twindragon by (— |, 0, 0)' so that its center is 0. □ 

Example 6.3. ( The non-fractal twindragon A ) 

A , the self-affine cube with 2 pieces, is the three-dimensional analogue of Example 
2.1. When considered as a rectangular parallelepiped with side lengths 1, \/2, ^4, it 
is the self-similar exception mentioned in Theorem 2.3, (ii). The similarity map g 
then is a 120° rotation composed with a homothety by the factor Of course, this 
tile has 26 neighbors: 6 faces, 12 edge neighbors and 8 point neighbors. 

Number of neighbors of twindragons. The twindragons form a nice little 
family with seven examples on which we can test our neighbor methods. However, 
the six fractal examples are truly self-affine, which makes their topology intricate. 
As a first overview, we compare the moduli of eigenvalues of M - the inverses of 
the contraction factors of the f) along different directions - and the number of 
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Figure 5. Twindragon £ 

neighbors. If the contraction factors differ much, there are lots of neighbors. We 
also determine the number of infinite boundary sets, by deleting from the neighbor 
graph the vertices representing finite boundary sets, as was explained in Example 
4.4 and Theorem 4.5. It turns out that all infinite boundary sets are uncountable. 
The number of faces will be determined in the sections below. 



Proposition 6.4. (First neighbor calculation for the seven twindragons) 



Twindragon 


A 


B 


c 


V 


£ 


JF 


g 


Parameter 


(0,0) 


(-1,1) 


(1,-1) 


(0,1) 


(2,-2) 


(1,0) 


(0,2) 


Complex eigenvalues ' modulus 


1.26 


1.29 


1.22 


1.15 


1.14 


1.09 


1.06 


Real eigenvalue 


1.26 


1.21 


1.35 


1.52 


1.54 


1.70 


1.77 


Number of neighbors 


26 


18 


20 


34 


34 


48 


76 


Uncountable boundary sets 


18 


U 


12 


14 


32 


48 


76 


Faces 


6 


14 


12 


14 


12 


16 


> 22 



For B , C and T> the number of infinite boundary sets is smaller than for the cube, 
so that we can expect them to be more or less ball-like. £ has medium complexity, 
and the last two examples have many infinite neighbor sets - they could be three- 
dimensional counterparts of Example 4.4. 

When we compare with the figures, we see that the twindragons possess a kind 
of fibre structure along the direction of the smaller eigenvalues of M, that is, the 
larger eigenvalues of M~ l . For B , the pieces drawn in Figure H] are stretched along 
the direction of the real eigenvalue of M~ l which is the largest contraction factor. 
So roughly the pieces look like 'cigars'. In all other examples, the larger contraction 
factor is given by the complex eigenvalues, and the pieces look like leaves or plates. 
For C in Figure HJ the leaves still have some thickness, while for £ illustrated in 
Figure [5j the plates are already very thin. For the other three twindragons, the 
plates look similar or still thinner. 
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a 


b 


c 


d 


e 


f 


g 


h 


i 


j 


1 





-1 








1 


1 


-1 


1 


2 





1 


1 





-1 


-1 








-2 


-1 











1 


1 


1 


1 


1 


1 






FIGURE 6. The neighbor graph for C contains all information about 
the topology of Figure [TJ The translation vectors, which we shall not 
need, indicate the position of the center of the neighbor. 



How to draw neighbor graphs of twindragons. Let us mention some other 
nice properties of our family and fix some notation for the neighbor graphs. First, 
all twindragons are connected. For two pieces this follows from Ti fl T 2 ^ (see 
Barnsley [13], 8.2.1). If this relation would not hold, T would be a Cantor set and 
would have no interior points. 

Next, the only possible differences kj — ki are k 2 — ki = (1,0,0)', k\ — k 2 = 
(— 1, 0, 0)', and ki — ki = 0. The corresponding labels i of the edges are 1 and 2 and i 
where the last case means two edges labelled with 1 and 2. We shall draw edges for 
the last case as double arrows, and we distinguish labels 1 and 2 by assigning a fat 
tip to the arrows with label 2. This way we save the labelling of edges. The label 
of the opposite edge is j = 3 — i for simple arrows, and j — i for double arrows. 

In view of the symmetry, we can simplify the neighbor graph further by drawing 
only one vertex from each pair of opposite vertices k, —k and introducing the con- 
vention that a wavy arrow ~» to k denotes an arrow to —k. This will also simplify 
the equation systems (jU]) for the boundary sets. It really matters whether one has 
to work with 20 or 40 vertices! Note that the opposite vertex of k is — k not only as 
a vector, but by Proposition 6.1 also as a boundary set: = —B^. 

For sake of brevity, we shall provide no details for the complicated twindragons 
J 7 and Q , but we give full arguments for the other four examples. 

To verify the condition of Theorem 5.1, we check in Figures El [TJ El and |9] that 
each vertex is reached either by a double arrow, or by two ordinary arrows with 
different tips indicating label 1 and 2, or by an ordinary and a wavy arrow with 
equal tips. In the last case the wavy arrow leads to —k, so we have to count the 
opposite edge. 
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FIGURE 7. Neighbor graph for the twindragon B . There are two 
point neighbors: i with address 2 and h with address 12. 

Example 6.5. We derive the equations (TJJ) of the boundary sets from Figure® We 
neglect all point neighbors, given by vertices below the thin line. 

B a = h(B b ) U f 2 (B h ) U f 2 (B c ), B h = f\{B A ) U f 2 (B d ), B d = f 2 (B { ), 

B c = MB,) U f 2 (B e ) U A(B f ), B c = f 2 (-B c ), B t = f 2 {B & ) . 
By substitution we reduce the system to two equations: 

2 2 

£a = h{B c ) U (J ^ 22 (i?a) , B c = f 12 (B a ) U |J f i2 (-B c ) . 

i,j=l i=l 

7. Identifying faces: the simple case 

Defining faces. The topological boundary of a set T will be denoted by 
dT — T \ int T. If T is a tile in a tiling, then <9T must be covered by the boundary 
of neighboring tiles: dT = [J k&v B k . 

Theorem 4.5 says that we can determine the cardinality of a boundary set B k by 
counting the number of paths starting from k in the neighbor graph G. It is much 
more complicated to determine the topology of B k . In three-dimensional self-affine 
tiles, we would like to distinguish boundary sets which are two-dimensional fractals, 
one- dimensional curves, Cantor sets or just finite sets. 

Let us say that the boundary set B k = T D (T + k) is a face of T if there is a point 
x G Bp, and an open neighborhood U of x with U fl dT C B k and U C T U (T + k). 
In that case, all boundary points of U H intT must belong to B k . Thus B k has 
topological dimension 2 since the boundary of each open set in R n has dimension 
n — 1. The following statement shows that we have six faces in Example 4.4. 

Proposition 7.1. Let k ^ be a vertex in G such that for each other vertex k' ^ 0, 
either k' or the opposite vertex —k' can be reached by a path from k. Then B k is a 
face ofT. 
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FIGURE 8. Neighbor graph for the twindragon T> . The neighbors 
i,r,p and k describe finite boundary sets with more than one point. 



Proof. Each tile must have faces. By fl9]), any path from k to k' implies the 
existence of a small copy of By in B^. Thus our assumption implies that B^ contains 
copies of faces, so it must be a face (cf. Proposition 8.1). □ 

For twindragons, there is always an edge with label 1 from the root to the first 
other vertex a = k 2 — k\ = (1,0,0)'. and this must be a face by Proposition 7.1. 
It turns out that for the twindragons B , C and T> , all infinite boundary sets fulfil 
the assumption of the proposition since from their vertices in G, there is a path to 
a. The neighbor graphs are shown in Figures [fU [7J and [HJ The point neighbors are 
separated from the faces by a thin line. They are found by Theorem 4.5, (iii). 

For B this can be seen in Figure [7J For C we have the loop at g with address 
2, and the 2-cycle {i, h} in Figure [6] which denotes a 4-cycle {i, — h, — i, h} in the 
complete neighbor graph, with address 2112. (The wavy arrow indicates that we go 
to the opposite sets, and there the labels 1 and 2 have to be interchanged. A second 
wavy arrow reestablishes the original labels.) In Figure El we have a similar 4-cycle 
{q, s,j,n} which denotes an 8-cycle of point neighbors in the complete neighbor 
graph with address 21221211. Also h and m are point neighbors, while i,r,p and k 
describe finite boundary sets with more than one point. We summarize: 

Proposition 7.2. (Faces of B , C andV ) 

The twindragon C has 12 faces and 8 point neighbors. B has 14 faces and 4 point 
neighbors. The twindragon T> has 14 faces, 12 point neighbors and 8 further finite 
boundary sets. 
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8. Different types of boundary sets 

The partial order for boundary sets. To compare different infinite boundary 
sets, we define an order on the vertices of G. For u, v G V we write u y v if there 
is a directed path from u to v. We say that the vertices are equivalent, and write 
u ~ v if either u = v or both u y v and v y u. For the boundary sets, u y v means 
that B u D fi(B v ) for some word i = i\...i q . 

Proposition 8.1. If u y v then dimB u > dimi?„ holds for the topological dimen- 
sion, the Hausdorff dimension and the box dimension. Moreover, if B v is a face, 
then also B u will be a face. 

Proof. The dimension concepts mentioned here (and also some others [17] ) have 
the following properties. 

(i) If A C B then dim A < dim B. 

(ii) If h is an affine map then dim h(B) = dim B. 

Together with the remark above this proves the first assertion. 
If B v = TC\T + v is a face, there is x G B v with an open neighborhood U G TUT + v. 
Now u y v means that B u D fi(B v ) for some word i where T + u D fi(T + v). So 
V — fi( x ) and its neighborhood fi{U) can be taken to show that B u is a face. □ 

Strong components of the neighbor graph. We can now classify the ver- 
tices of G in a similar way as states of a Markov chain. Clearly ~ is an equivalence 
relation. The topological and Hausdorff dimension, as well as the property of de- 
scribing a face, are invariant with respect to ~ . Let V be the set of equivalence 
classes with respect to ~ which are called strong components of G. Then we get a 
new graph G by drawing an edge from (u, v) if there is u G u and v G v such that 
(u, v) is an edge in G. 

The graph G compresses the structure of G. If two vertices u,v G G belong to 
the same class, this means that B u contains a copy of B v and conversely, so both 
boundary sets have essentially the same structure. The root G V would also be 
the root of G but it has no meaning and is omitted. 

The relation >- is a partial order on V. If u y v then all boundary sets B u contain 
copies of all boundary sets B v , and dim£? u > dim.Ru for all choices u G u,v G v. 
Thus the largest classes of boundary sets are the faces, then we have one-dimensional 
sets, Cantor sets, and as terminal classes countable and finite sets. Boundary sets 
from classes u, v with u y v can also have equal dimensions, however. 

Remark 8.2. Consider a class u EV, as a subset ofV. There are three cases. 

(i) u contains only one vertex. Then the boundary set B u has the same dimen- 
sion as the successor vertex v (in case of several successors, the union of the 
successor boundary sets). Namely B u = fi(B v ). 

(ii) u is a cycle in G, without loops and diagonals. In that case the boundary sets 
B u with u G u are countable unions of copies of the successor boundary sets. 

(iii) u contains two cycles. Only in this case u produces its own Cantor structure, 
and the B u can have larger dimension than the successor sets. 
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In Figures El d and El we have one class of faces, a number of singleton classes 
and a terminal cycle which is a class of two, one or four points. In Figure Owe have 
a class of six faces, and a class of eight infinite neighbors h..r on the right, as well as 
three singleton classes. The vertex s describes a point neighbor, the vertices g and 
p correspond to infinite boundary sets with equal dimensions as h..r, by the above 
remark, (i). We now have to decide whether g,p and h..r are faces. 

9. Identifying faces: the general case 

An algorithmic approach to find the faces. When we consider one equiv- 
alence class of the neighbor graph which contains a face, then it consists of faces 
only. 

Problem 9.1. Must all faces belong to one equivalence class of G ? 

This problem also arises when we work with the contact matrix. In [2\ , irreducibil- 
ity was just assumed. We present methods to decide the problem algorithmically. 
We start with a simple topological observation. 

Proposition 9.2. The faces cover the boundary of a self-affine lattice tile T. 

Proof. Let i be a boundary point of T and U an open neighborhood of x. Let 
E = U fl dT and let T k , k = l,...,q denote the neighbor tiles of T. The sets T fl T k 
cover E, and by the Baire category theorem one of these sets, say TDTj, must 
contain an open subset V of E. Each y G V has a neighborhood W in MJ 1 which is 
contained in T U T\. So T fl T\ is a face. Since this holds for any U, the point x is 
an accumulation point of the faces. Since we have only finitely many faces, x must 
belong to one of them. □ 

Theorem 9.3. (The neighbor graph decides which boundary sets are faces) 

Let us assume that all neighbors of the tile T are compatible. Then a boundary set 

B u is a face if it is not contained in the union of all other boundary sets. 

For the neighbor graph this means that there is a word i which belongs to the language 

L u of labels of paths starting in u, but not to L v for any other v ^ u. 

Proof. Consider the union of the tile T and all its neighbor tiles. If B u is not a 
face, then by the proposition it is contained in the union of the remaining boundary 
sets. If B u is a face, there is x G B u and a neighborhood U cTUT + uofx. Then 
U cannot intersect any other neighbor T + v since neighbors must have no common 
interior points. The first part of the theorem is proved. 

Now x has an address: x = 7r(fi) where 7r is the projection from the symbol space 
I°° to T. Since 7r is continuous [HI [17] , there is a prefix i of x such that all addresses 
t which begin with i will fulfil 7r(t) G U. By Theorem 4.5 (iii), this shows that for a 
face B u there is a word i which is contained in L u and in no other L v . 
Conversely, suppose B u is not a face and i G L u . There is an address fi which begins 
with i such that 7r(B) = x is not on the intersection of different pieces fj(T). (It is 
known that the intersections have Lebesgue measure [TJ, and since the normalized 
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FIGURE 9. Neighbor graph of twindragon £ . The strong component 
of a..f represents faces, s is a point neighbor. Vertices h..r form a 
strong component representing infinite boundary sets, g and p are 
components with boundary sets of the same dimension as h..r. 

Lebesgue measure on T can be considered as the image measure of the product mea- 
sure {—,...,—} on J°°, no cylinder can be mapped completely into the intersections.) 
Since x is also in a face B v and has only one address, L v contains the sequence B 
starting with i. □ 

Examples of address calculations for boundary sets. We shall apply our 
theorem to Figure |9j Address calculations are boring and best left to the computer, 
but they can provide valuable information about the intersection of pieces. We start 
with simple examples. 

Example 9.4. ( Point neighbors in C and B ) 

Which faces contain the point neighbors of C and El In Figured we find 

B s = tt(2) = B a n B h n B c n B d , B h = 7r(22TT) C B c n B d n B c and 

Bi = 7r(2Tl2) = fi_ c nB c nB f , B s = 7r(2l22T) C5 a n B_ c 

so g,h,i seem proper vertices but j seems to be on an edge. On the other hand, both 
point neighbors of Figure\7\ coincide with the intersection of only two faces: 

Bi = tt(2) = B c n B c , B h = tt(12) = B a n B c . 

Already these calculations show the strange polyhedral structure of twindragons. 
For Figure [9] the boundary sets look much more complicated than for ordinary 
polyhedra where an edge is covered by two faces. The proof of the next proposition 
will also reveal hidden symmetry behind the apparently orderless structure of G. 

Proposition 9.5. ( Boundary structure of £ ) 

The vertices g, p, h..r in the neighbor graph of the twindragon £ do not represent 
faces. The boundary set i is contained in the union of the four faces e, -f, c, -a and 
not in a union of three of them. 
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Proof. We determine the language Lj from the three elementary cycles which 
lead in G from i to either i or -i (cf. section 7 and Example 9.4). 

4 

Li = *2111L i U 12L_; U * * *1122L_i = C X U U C 2 U U (J D i L_ i (11) 

3=0 

where * is considered as a wildcard for both 1 and 2, and Ci = 12111, C2 = 
22111, D = 12, D x = 1 * 11122, D 2 = l* 21122, £> 3 = 2 * 11122, £> 4 = 2 * 21122. We 
now consider the set of faces J = {e, — f , c, —a}. We will show 

U C Lj = L C U L_ f U L c U L_ a . (12) 

First we study the action of C\ and C 2 on the vertices of J. In the following table 
+ means there is path from x to y with label C\, and * means the same for C 2 . 



x y 


e 


-f 


c 


-a 


e 


* 




* 




-f 




+ 




+ 


c 




* 




* 


-a 


+ 




+ 





Claim 1: any word C G {C\,C 2 }* is in Lj, and it can be realized with any 
prescribed terminal point in J. 

This claim can be proved by induction on the number t of terms in C = (7™ 1 ...C'™ t , 
using the table. The starting point always depends on the word. For instance, it 
must be -f for C = Cf. 

Now we study how the Dj, j = 0, ..,4 lead from J to — J. In G it can be seen that 
D leads from e to -e, -f to f, c to a, and -a to -c, respectively. Moreover, D\ leads 
from -a to all vertices of —J, and from no other vertex of J to any vertex of —J. 
Similarly, D 2 leads exclusively from -f to all of — J, and D 3 leads exclusively from e 
to all of —J, and D4 leads exclusively from c to all of — J. From this observation it 
follows that all vertices of J are needed to obtain L x C Lj. Moreover, it is clear that 

Claim 2: any word Dj,j = 0, ..,4 can be realized by a path from J to — J with 
any prescribed terminal point in — J, provided we can select the starting point. 

Applying Claims 1 and 2 alternatingly, to J and its opposite set — J, we complete 
the proof of f|T2|) . The rest follows from Theorem 9.3. □ 

10. Dimensions and eigenvalues 

The modified Hausdorff dimension. In this section we give an alternative 
proof for the fact that £ has only 12 faces. It is based on recent work of He and 
Lau [23] and of Akiyama and Loridant [2]. Let H denote the adjacency matrix of 
the graph G \ {0}. For two non-zero vertices u,v the entry h uv denotes the number 
of edges leading from u to v. It is also possible to take the adjacency matrix of our 
simplified graphs, where only one vertex from each pair {k, —k] is taken into account 
and h uv counts the edges from u to {v , — v}. Since H is a non-negative matrix it has 
a real eigenvalue of maximum modulus which we call A. 
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When the matrix M defining the affine maps fj is conjugate to a similarity matrix 
with expanding factor R (see Proposition 2.2) it was shown by Mauldin and Williams 
[34] that the Hausdorff dimension of the sets from (Q fulfils 

dim5 fc <-^A, (13) 
log R 

and equality holds if an open set condition is satisfied, see the discussion at the end 
of this section. 

For self-affine sets, however, calculation of dimensions is much more complicated 
then for self-similar ones. However, for our case where all fj are defined by one ma- 
trix M, He and Lau [23] have defined a pseudo-norm w on M. n which turns the fj into 
similarity mappings, with factor r = 1/ y det M. Moreover, w fulfils a weak type of 
triangle inequality which is sufficient to develop Hausdorff measure and dimension 
as usual. And w generates the Euclidean topology even though it drastically distorts 
the geometry of MJ 1 . The value of the modified Hausdorff dimension dim™ B u was 
used in [23J to give estimates of the real Hausdorff dimension in terms of eigenval- 
ues of M. Akiyama and Loridant [2] applied the modified Hausdorff dimension to 
boundary sets of tiles. Here we shall compare it with the topological dimension. 

Theorem 10.1. (Upper estimate of topological dimension of boundary sets) 
Suppose that T is a self-affine tile with m pieces in M. n , and all neighbors are com- 
patible. Let W denote a strong component of the graph G, and let Xw denote the 
largest eigenvalue of the adjacency matrix Hw of W as a subgraph of G. If for some 
integer q < n, 

n log Xw 

-j < q 

logm 

then the topological dimension of all boundary sets B v with v G W is strictly smaller 
than q. 

Proof. The matrix Hw is irreducible, and W is non-periodic (i.e. the g.c.d. of the 
lengths of cycles from a point v G W to itself is 1) because of the compatibility of 
neighbors. So Xw is the unique largest eigenvalue. The Hausdorff dimension of the 
sets B v with v G W with respect to the He-Lau pseudo-norm is < = n 

since R = yfm [21 [23] . The pseudo-norm generates the Euclidean topology, and the 
topological dimension is always smaller than the Hausdorff dimension. □ 

Example 10.2. ( Boundary sets of the twindragon £ ) 

For S we determine the dimension of the boundary sets B v with v G {g, p,h..r}. We 
obtain Xw ~ 1.554. A small polynomial for Xw can be obtained from equation ( fTIj) 
for B = B x which expresses B as a self-affine set 

2 2 

B = \Jf t2U i(B)Uf 12 (-B)U [j f ijkll 22(-B) . 

i=l i,j,k=l 

The characteristic equation is 1 = 2z 5 + z 2 + 8z 7 , which for X = 1/z becomes 
X 7 - X 5 -2X 2 -8 = 0. 
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The modified Hausdorff dimension is log Xw / log v^2 ~ 1.908. Thus the topological 
dimension of the B v is < 1. So the B v cannot be faces. A calculation for the faces 
a..f gives modified dimension 2.67. 

For the complicated twindragons J 7 and Q we also computed the strong compo- 
nents and their modified dimensions. For Q , the above theorem does not apply 
since beside the irreducible component of 22 faces with modified dimension 2.87, 
there is another component of 20 boundary sets, plus 2 singleton components, all 
with modified dimension 2.13. All the remaining 32 boundary sets have modified 
dimension 1.52. This boundary seems really complicated! The other example is less 
intricate. We omit the details of the following example. 

Example 10.3. ( Boundary sets of the twindragon J 7 ) 

For J 7 , the neighbor graph G \ {0} has an irreducible component V\ of 16 faces 
with modified dimension 2.78, and the other 32 neighbors are not faces. 10 of them 
form a cycle without diagonals as second component V2, and there are two singleton 
components between V\ and Vj>. Since the cycle Vi consists of six double arrows and 
4 simple arrows, A2 = 2 3 / 5 and 3 * log A2/ log 2 = 1.8. Due to the theorem, the 
topological dimension is at most 1. 

Moreover, there is another 10-cycle V3 with only two double arrows and A3 = 2 1 / 5 , 
and 10 singleton components between V% and V3. So the 20 corresponding boundary 
sets have modified dimension 3 * log A3/ log 2 = 0.6 and by Theorem 10.1 must be 
Cantor sets. 

When one considers the labels along the cycles V2 and V3, one obtains the languages 

L 2 = *1 * *1 * 2 * *2 and L 3 = 11 * 1122 * 22 

and their cyclic permutations as address sets of the corresponding boundary sets. 
From this it is easy to conclude that each Cantor set associated to a vertex in V3 is 
the intersection of two respective boundary sets in V^. 

We did not use the open set condition of the boundary equations (jHJ), and this 
condition need not always be fulfilled. In the cube (Example 6.3), for instance, 
some of the equations © contain equal terms since an edge is counted two times. 
However, the open set condition will be true if we remove such double sets and if 
the neighbors are all compatible. In particular it is true when we restrict ourselves 
to faces [2]. When the graph is a cycle, as V2 and V3 in the last example, the open 
set condition also holds. 

11. Polyhedral structure of tiles 

Homeomorphism with a polyhedron. We have seen that the structure of 
boundary sets can be quite complicated. We are interested in the presence of simple 
structure, however. Let us say that Tel 3 has the structure of polyhedron if T is 
homeomorphic to a ball in M 3 , the faces are homeomorphic to disks in IR 2 , and their 
non-empty intersections are either singletons or homeomorphic to an interval - in 
such a way that T together with its boundary structure becomes homeomorphic to 
a polyhedron in Euclidean M 3 . 
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This condition is stronger than just requiring that T is homeomorphic to a ball. 
However, this is the structure which crystallographers expect of a tiling: "We define 
a tiling as a periodic subdivision of space into bounded, connected regions without 
holes, which we call tiles. If two tiles meet along a surface, we call the surface a 
face. If three or more faces meet along a curve, we call the curve an edge. Finally, 
if at least three edges meet at a point, we call that point a vertex" |18j . 

Definition of edges and vertices. We follow this quotation and define a 
polyhedral structure on every self-affine tile. We remove from G all vertices which 
do not represent faces. Then we define the graph G 2 of edges - that is, of intersection 
sets Bk D By. For these, it is not so difficult to check whether they are homeomor- 
phic to intervals [TO]- If necessary, we can also determine the graph G 3 of vertices 
of the tiling. Euler's polyhedra formula can then be used to check whether T is 
homeomorphic to an ordinary polyhedron. The method is quite fast in showing that 
certain tiles are not ball-like. 

A definition similar to the following was given by Thuswaldner and Scheicher [36] 
but it was not applied to examples. 

Definition 11.1. Let T C M™ be a self-affine lattice tile with neighbor graph G = 
(V, E = UZLi Ei)- F° r ^ = 2, 3 the graph of £ -intersections of boundary sets ofT is 

G l = (V\ E e ) with V 1 = {K = {k x , h} | all k { different, B kl n ... n B k( ^ 0} 

rn 

and E £ = {jE e i with (14) 

i=X 

E- = {(K, K', i) | there is a 1-1-map : K -> K' with (k, (j)(k),i) G E { for k 6 K} 

The idea is that a neighbor intersection (T + k) H (T + k') has a piece which is 
at the boundary of Tj if and only if both neighbor sets have such a piece, and the 
pieces do intersect. The algorithmic check for non-empty intersection is the same as 
for G. We start with all ^-subsets of V, and step by step we remove those which have 
no outgoing vertices to the remaining vertex set. Since for £ = 2, 3 we must have 
non-empty intersections, we end with a graph where each set has outgoing vertices, 
and hence 'infinite paths' of edges (counting repeated use of cycles). The properties 
of G l are similar to those of G. 

Proposition 11.2. G l provides a graph- directed fractal structure on the £-intersec- 
tions of neighbors. For each vertex K e V , the labels of infinite paths in G e 
with starting point K are the addresses of the corresponding boundary intersection 
f) keK B k . The cardinality of this intersection can be determined as in Theorem 4-5. 

Example 11.3. We continue Example 4-4 by determining G 2 for the six faces of 
the Sierpinski tile, see Figure^ Since the outgoing edges from —a, — b, — c all start 
with 4, 2, and 3, respectively, the corresponding boundary sets have pairwise empty 
intersection. The intersection of a boundary set with its opposite is also empty. We 
get edges 

{a, — b} A {a, — c} A {a, — b} 



26 



CHRISTOPH BANDT 



so S a ni?_b is a singleton with address 2323... = 23 and B a (~)B^ c has address 32. The 
vertices a, b, c all have a loop with label 1, so the corresponding three boundary sets 
meet in the point with address 1. Altogether, G 2 consists of three cycles of length 2, 
and three isolated loops. G 3 has only one vertex with a loop, representing the center 
point of the tile. 

12. Polyhedral structure of twindragons 

Intersections between antipodal surfaces. The first information which the 
face intersection graph G 2 provides is a list of non-empty intersections of pairs of 
faces. As in Theorem 4.5, it is easy to decide whether such an intersection consists 
just of one point. 

Proposition 12.1. Two faces B u , B v with u,v &V intersect each other if and only 
if {u, v} is a vertex in G 2 . Moreover, if there is only one infinite directed path in G 2 
starting in the vertex {u, v}, then B U DB V is one point, and the address of this point 
can be read from the labels of the path. 

For our twindragons we check whether there is a face Bk which intersects its 
opposite face B^. If such intersections exist, it is unlikely that our tile has the 
structure of a polyhedron. Actually, this happens for four of our twindragons and 
provides an argument to show that T> and £ are not homeomorphic to a ball. 

A set A C W 1 is called simply connected if it is connected and its homotopy group 
is trivial. The latter means that A has "no holes" - any closed curve C C A can be 
contracted within A to a point. 

Theorem 12.2. (Tiles for which the interior is not simply connected) 

Let T be a connected self-affine lattice tile which has a symmetry center c which is 

not an interior point of T. Then the interior of T is not simply connected. 

In particular, the interior of V and S is not simply connected. These twindragons 

are therefore not homeomorphic to a ball. 

Proof. We can assume c = and thus T = — T by changing the maps as in 
Proposition 6.1. If is not an interior point of T, there is some tile T + k in the 
tiling which contains 0, and by symmetry also — (T + k) = T — k does not intersect 
intT and contains 0. Take a big sphere S around which contains the three tiles 
in its interior. Let z G S be of the form z — t • k, t > 0. There is a line segment a 
from z to the point y G T + k which maximizes the scalar product < y,k > . Since 
T is connected, it is arcwise connected, and there is an arc (5 C T + k from y to 0. 
The union e of the arcs a, (3, —(3, —a is an arc in S from z to —z which does not 
intersect int T. 

Now assume int T is connected and take a point x G int T. The opposite point 
—x is also in int T. So there is a connecting arc from x to —x. Covering this arc by 
finitely many open balls inside int T, we see that we can replace it by a union 7 of 
finitely many line segments, and we can choose them so that no two segments are 
parallel. Then the union 7 U —7 is a closed curve from x through —x to x with a 
finite number of self-intersections. If x\ is the first self-intersection point, we find a 



COMBINATORIAL TOPOLOGY OF THREE-DIMENSIONAL SELF-AFFINE TILES 



27 



closed arc from x\ through —x\ to x\ with two self-intersection points less. After a 
finite number of steps we get a simple closed curve 5 in int T from some x n through 
— x n to x n . By construction, 5 surrounds the arc e and thus cannot be contracted 
within int T to a single point. So int T is not simply connected. The second assertion 
will be proved in the example below. □ 

Proposition 12.3. If an infinite path in the neighbor graph G of a twindragon starts 
in the root and contains no double arrow, it defines an address of the center ofT. 

Proof. For a path of single arrows which addresses a point x, the opposite path 
is obtained by just interchanging labels 1 and 2. In the proof of Proposition 6.1, it 
was shown that the address of the opposite path corresponds to the point —x. For 
paths starting in the root of G, Theorem 4.5 (i) says that both addresses belong to 
the same point, x = —x implies x — 0. □ 

Example 12.4. ( Faces intersecting their opposite face in T> and £ ) 
For T> , Figure [21 shows only one path from the root without double arrows, given 
by the cycle acfd which describes an 8-cycle in the complete neighbor graph. The 
associated address 122211112 corresponds to the center 0. As can be directly seen 
in Figure 0, this address is also obtained from a path starting in vertex c , going 
through the cycle cefg and the opposite vertices. Thus the center belongs to the 
boundary set B c , hence not to the interior ofT. 

For £ , two root addresses of the center are visible in Figured 1212 from abs, 
and 1221221 from abim. The second address is also obtained when we start at -e 
and go -ebcef-e-fef... Thus G B_ e , the center is not in intT. 

The addresses were found from a calculation of G 2 although G is sufficient to 
check them. We determined G 2 also for the more complicated twindragons and 
found that for J 7 there are three faces B^ which intersect their opposite faces -B_&, 
not in a point, but in a Cantor set with dim™ B^DB^ = 0.6. It is not clear whether 
any of the Bk contains c, so the above argument fails, but it is obvious that T has 
not the structure of a polyhedron. In an ordinary polyhedron, given topologically as 
a planar map on the sphere, at most one pair of opposite faces can intersect. For Q 
the situation is similar: we have even five pairs of opposite faces Bk which intersect 
each other, not in a point, but in an uncountable set. 

The polyhedral structure of B and C . Only the two twindragons shown in 
Figures [1] and H] can still have the structure of a polyhedron. We shall prove that 
this is not quite the case but we expect that 

Conjecture 12.5. The twindragons B and C are homeomorphic to a ball. 

It will not be possible to settle this question here, but we shall establish the 
polyhedral structure which leads us to the conjecture. We start with the list of 
non-empty intersections of faces which was established by computing the graph G 2 . 
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Figure 10. Polyhedral structure of C (left) and B (right) 

Example 12.6. ( Polyhedral structure of C ) 

In C , the following faces have uncountable intersection: 

a with b,f,-c,-e b with a,c,d,-e c with b,d,-a,-f 

d with b,c,e,f e with d,f,-a,-b f with a,d,e,-c 

Moreover, there are one-point intersections 

a,-a with d,-d, b,-b with f,-f and c,-c with e,-e 

and two-point intersections a with c and -a with-c. The two points coincide, however, 
since the corresponding addresses are identified by Theorem 4-5 (i). 

Proof. The proof is a simple but lengthy calculation. We sketch some facts which 
can be seen directly from inspection of G in Figure |6j The address 2 belongs to 
L a , L b , L d , Lf and it is clear that it is the only address in L a nL d as well as in L b flL f . 
For the opposite faces we have address 1. Moreover, L c D L e = 2211 = L d fl L_ a 
and L c fl L_ e = 1221 = Lf fl L_ b . We can consider the languages restricted to the 
vertices of faces, due to Proposition 9.2. Nevertheless, we see that these addresses 
also appear for the point neighbors g, h and -j, respectively, so they should really 
represent corner points in the tiling by T. 

It is obvious from G that Lf fl L_f = Lf fl L_ c = 0. From this we conclude 

L b n L c = 2(L d n L_ c ) = 22(L f n L_ f ) U 22(L f n L_ c ) = 0. 

Now we can determine 

L a n L c = *(L b n L e ) u 2(L C n L e ) u i(L b n L f ) = 222TT u 12 , 

and these two addresses are equivalent since they label opposite paths from the root 
ofG. ' ' □ 

Proposition 12.7. ( Equivalent convex polyhedron for C ) 

The polyhedron on the left of Figure UU realizes exactly all intersections of the faces 
ofC , except for B & C\B C and L?_ a nL?_ c . It is an octahedron truncated at four vertices, 
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in the middle of the corresponding edges, and its faces are rhombi. 

Nevertheless, C has not the structure of an ordinary polyhedron because the interior 

of the face is not connected and so is not homeomorphic to a disk. 

Proof. The first part is done by checking all pairs of faces. For the second part 
becomes clear when we prove that B & fl B c lies on the surface of T, in B^. However, 
it can be seen in G that 12 is an address in Lb. 

Here is another argument. From G we have the equation B h = f\(B d ) U f2(B d ). 
On the other hand fi(T) n f 2 (T) = /i(B a ) = / 2 (S-»). We have seen that B a and 
£>_ a intersect B d in a single point, and this also holds for their images under fi or 
$2- Thus /i(.Bd) and f2(B d ) have only one common point. So £>b consists of two 
isometric pieces which intersect in a singleton, and cannot be homeomorphic to a 
disk. □ 

It should be mentioned that the polyhedron in Figure [TU] describes C only topo- 
logically - the metric structure is completely different. We have B{ = /^(-Ba) and 
Bd = f2{B{) so these faces are of different size. The intersection of the faces 
£ a , B h , B d , B { is not a corner as in FigureEUJ the point tt(2) = (|, 0, -)' = 2±^ = ^ 
is on the middle between the centers of opposite adjoining faces (see the table at 
Figure [6]). A similar calculation shows that B & C\B C is in fact the center point of B^. 
Now let us study B in the same way as C . 

Example 12.8. ( Polyhedral structure of B ) 

The following faces have uncountable intersection in B : 

a with c,d,e,g,-b,-f b with c,f,-a,-d c with a,b,e,f,g,-d 

d with a,e,f,-b,-c,-g e with a,c,d,f f with b,c,d,e,-a,-g g with a,c,-d,-f. 

Moreover, there are one-point intersections 

B c n -B-c = 7r(12), B c PI -B_ a = 7r(221), and their opposites. 

Proposition 12.9. ( Equivalent convex polyhedron for B ) 

The polyhedron on the right of FigurelTdi realizes exactly all intersections of the faces 
of B , except for the point neighbors. It is an octahedron truncated at all vertices, at 
one third of the corresponding edges. The faces are regular hexagons and squares. 
Nevertheless, B has not the structure of an ordinary polyhedron because the interior 
of the faces B^ and B d is not connected and so these faces are not homeomorphic 
to a disk. 

This is proved similarly as for C . Note that faces of the truncated octahedron 
never meet in a single point. Using Figure [7] we see that B c fl 5_ a is on B^ and 
B c n B_ c on B_ d . 

13. Connectedness of the interior 

In this section, we go a small step towards proving our conjecture, by proving 
that the interior of C is connected. 

Proposition 13.1. Suppose T = \JuLxfj(T) is a self-affine tile and there is a 
connected set E C intT such that E fl fj{E) ^ for j £ I = {1, ...,m}. Then the 
interior of T is connected. 
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Proof. D\ = UjeJ fAE) * s a connected subset of int T. By induction we show that 
D n = Ukui<n fw{E) is connected, where \w\ denotes the length of the word w. Then 
D = U^°=i D n = \J we j* fw{E) is connected, and this set is dense in intT. Hence 
int T is connected. □ 

Proposition 13.2. intC is connected. 

Proof. A piece T w belongs to int T if it is surrounded in T by all possible 
neighbors. T w has all possible neighbors within T if a path labelled with a suffix 
of w leads from the root to each vertex of G (cf. proof of Theorem 5.1). When we 
require this for faces only, w = 121212 and v = 2121121 fulfil this condition, as one 
can check with Figure El We put 

E = T' w U T' v U % U % , w = 212121, v = 1212212 

where T' w means T w minus its one-point boundary sets. Then E C int T. Since E 
contains 7r(12) and 7r(21), the set E intersects fi{E) and f2{E). To apply Proposition 
13.1, it only remains to verify that E is connected. Since twindragons are connected 
(cf. section 7), each T' u is connected. Moreover, there are pairs of opposite paths 
from the root in G with labels vol and v , hence also with w2 and v, and also with v 
and v, which all end in a vertex of G corresponding to a face. Thus T wl and T v , T% 
and Tu,2, as well as T v and Xg have a face in common. So E is connected, and the 
proposition applies. □ 

w and v were obtained in a straightforward way. w is the shortest word for which 
T w has all possible face neighbors (we used a computer search), v is the shortest 
word with all face neighbors which contains the center point by Proposition 12.3 
(this is easy: the address of is 1212). In fact, is an interior point of T v U T$. For 
B we also determined such v. The address of is 1122, which gives v = 1122122122. 
This indicates that the interior of B is more fragmented than that of C and, to let 
the above method work, we would need a longer chain of pieces. 

14. Remarks on algorithms 

Use of computer. We have shown that the topology of self-affine tiles can be 
studied by rather simple methods: finite automata and regular languages. To answer 
concrete questions, however, we have to perform quite a number of elementary logical 
operations. Computer assistance seems necessary and convenient to study more 
complicated examples. 

Although all results in this paper have been checked by hand, interactive computer 
work was essential to obtain them. Let us briefly document the main algorithms we 
have used. 

Construction of the neighbor graph. To construct G, three steps are 
performed, as indicated in Example 4.4. 

(i) Lower and upper bounds l q , b q for {x q \ x G T} are derived for each coordinate 
q = 1,2,3. They need not be sharp, but should be taken rather generously 
(we work with integers anyway). If b q is too large, the computer will work 
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a little longer - if it is too small, only part of the neighbor graph will be 
determined. For examples like B and C , 100000 points of an IFS algorithm 
[TTj are completely sufficient to find l q and b q . For J 7 and Q there are rare 
outliers on the thin fibres, and an exact estimate is needed. 

(ii) Starting with k = 0, a list of new vectors is calculated by the recursion 

k' = Mk + kj - h 

where i,j take all values in {1, ..,m}. The index of k is listed together with 
the labels i, j at k' as the predecessor of k', and the number of successors of 
k is updated at fc. Of course, k' is only processed if it is within the bounds 
l q and b q . Moreover, it must always be checked whether k' has already been 
listed before. In that case, the new predecessor is added at the old place. 
Since the number of possible k is finite, the whole list will be completed after 
finite time, and describes a graph containing G. 

(iii) All k without successors are removed from the list, and for each of their 
predecessors, the number of successors is updated. This step is repeated 
until each k has at least one successor (or, in case of a Cantor set, G is 
empty) . 

When the resulting G has predecessors of the root, the open set condition is not 
fulfilled. This cannot happen if the kj form a complete residue system for M. 

Other operations with graphs. For constructing the intersection graph G 2 
from G, we have to consider G x G, identify {u, v) with (v, u) and define the edges 
according to Definition 11.1. Then we apply step (iii) as above, to exclude empty 
intersections. 

To identify point neighbors in G, we use the adjacency matrix H defined in section 
10. If q is the number of non-zero vertices in G, then k is a point neighbor if the sum 
of the row associated with k in H, H 2 , H q always equals 1. When we repeatedly 
apply a procedure like (iii) to the point neighbors, we also remove the finite boundary 
sets. One-point intersection sets in G 2 are found in the same way as point neighbors 
inC7. 

The partial order of boundary sets (section 8) is also determined with H. Let the 
matrix N be the maximum of H + H 2 + ... + H q and 1, taken in each cell. Then 
n uv = 1 if and only if u >- v. In the row of u there are the elements which can be 
reached from u, in the column of u there are the v from which u can be reached. The 
strong components of G can be found by properly ordering the rows and columns 
ofiV. 

To deal with languages L u , we introduce adjacency matrices Hi for the edges with 
label i. Then the matrix H1H2H2, for instance, tell us between which vertices there 
is a path labelled 122. This fact was used to find the proof of Proposition 9.5. 

There are many similar tools waiting to be developed and tested to help reveal 
the geometric structure of fractal tiles. 
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